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ABSTRACT 


This thesis presents the results obtained in computer simulations 
of the commutated network filter. The commutated network filter is also 
used as a bandpass filter in servosystems. A linear approximation to 


the transfer function is found to be usable for design purposes. 
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I. INTRODUCTION 


A commutated network is an electrical network, either active or 
passive, consisting of resistors and reactive elements which are 
switched or commutated in a periodic manner. In most practical cases 
the reactive elements are capacitors. There are several configurations 
of a commutated network used in the field of electrical engineering. 
The configuration used here is that of commutated capacitors used in 
the feed-back path around an operational amplifier. Thus by applying 
commutation action to the capacitors in a lag network, the low-pass 
frequency response of Pia network can be transformed into that of a 
highly selective bandpass filter. The center frequency of the bandpass 


characteristic will be the same as the commutation frequency. 


Ea 


II. COMMUTATED NETWORK CONFIGURATION AND 
COMPUTER SIMULATIONS OF THE NETWORK 


The fundamental element of commutated networks is the capacitor 


Shown diner iO6 2 1) Cc ‘ 
>; ae 
a t 


Figure 2.1 Commutated capacitor 
The commutating function, P, is a unit square wave of the same 
frequency as the input signal, e.. The following equations can be 


WrLttene for yiloeee . ik 


l e 
e = at dt 
Cc Cc Cc 
i= Pei, 
Cc i 
ae SS oe 
Cc 


e = oa hae dt 
Cu i 
e = fi, at 
Cc i 
since P-P = 1 both when P is in the low and in the high position. 


Figure 2.2 shows how the capacitor is placed in the coupled RC 


Commutated network. T C ) 
“ene 





Figure 2.2 Coupled commutated network 
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(2-1) 
(2-2) 
(2-3) 
(2-4) 


(2-5) 





Another variation of the commutated network is the uncoupled RC network 


sliéwn in fig. 2.3. 


Loe 





Figure 2.3  Uncoupled commutated network 
Reference 1 proves that both variations of the commutated network 


can be represented by the circuit shown in fig. 2.4. 





Figure 2.4 Equivalent circuit of commutated networks 
The number of capacitors in the commutated network can be varied. 
For the following, four capacitors are used in the coupled version of 


the network. Figure 2.5 shows the circuit used in the analysis. 





Figure 2.5 Four-capacitor commutated network 
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The equivalent network for the four-capacitor coupled commutated network 


is shown in fig. 2.6. 





Figure 2.6 Equivalent circuit for figure 2.5 


A computer program was written - in IBM's Digital Simulation 


Language - to simulate the network of fig. 2.5. Figure 2.7 shows the 


program used. 


F = 1, 

X = SIN(F*2.*3,1415*TIME) 
ZONE = IMPUIS (.0,1.00000) 
ZTWO = IMPULS (.1250,1.00000) 
ZTHREE = IMPULS (.2500,1.00000) 
ZFOUR = IMPULS (.3750,1.00000) 
FONE = PULSE(ZONE, .499) 

FTWO = PULSE(ZTWO, .499) 
FTHREE = PULSE(ZTHREE, .499) 
FFOUR = PULSE(ZFOUR, .499) 
PONE = 2.*FONE -1. 

PTWO = 2.*FTWO -1l. 

PTHREE = 2.%*FTHREE -1. 

PFOUR = 2.*FFOUR -1. 

XONE = -X - .95%Y 

WONE = PONE*XONE 

WIWO = PTWO*XONE 

WIHREE = PTHREE*XONE 

WFOUR = PFOUR*XONE 

OONE = INTGRL(.0,WONE) 

OTWO = INTGRL(.0,WTWO) 

OTHREE = INTGRL(.0,WHTREE) 
OFOUR = INTGRL(.0,WFOUR) 
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YONE = PONE*(-OONE) 
YTWO = PTWO*(-OTWO) 
YTHREE = PTHREE*(-OTHREE) 
YFOUR = PFOUR*(-OFOUR) 
Y = -( YONE + YTWO + YTHREE + YFOUR) 
Figure 2.7 Computer program for simulation 
The optimal value of the feed-back gain, K, was used. Reference 2 
proves that the optimal value is 0.95 (1/(1+1/7 +1/9 +1/15 +1/17 + ...). 
The modulating square wave, p, was generated in the following way: 
an impluse, Zi» was generated at the beginning of every period. This 


impulse triggered a square-wave generator and the result, Fos is shown 


1 Nope ee tO - 


F 
4 > 2.8 
ne T+t. > +t 


Figure 2.8 Square wave, F. 
The square wave, Fs was multiplied by a factor of two and the DC-level 


was shifted minus one unit. The result is P. as showm in@fie. 2.9. 


Yves) 





Figure 2.9 Square wave, P 


The different square waves are delayed with respect to one another. P, 
is delayed T/8 with respect to Pi. P. is delayed T/8 with respect to Po» 
and P, is delayed T/8 with respect to P,. ty was chosen equal zero. 


Figure 2.10 shows the square waves as used in the program. 


1. 


The following equations can be written for the commutating waveforms. 


Py is a unit square wave of period T, starting at t=0 thus f = 1/T Hz 


and w = 2n/T. 
n 





Figure 2.10 Commutating square waves 
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The objective of the RC-commutated network is to construct a 


Ne 


signal of same frequency but of opposite phase as the input signal. 
If this is the case the sum of the input signal and the output signal 
has zero component of the input frequency. 


Figure 2.11 shows how this can be implemented. 


y 


X COMMUTATED 


NETWORK 








OO 








Figure 2.11 Filter configuration 

Here OUT is the sum of X and Y. As a result of the square waves being 
applied to the capacitors, harmonics of the frequency of the input 
Signal are present in the signal OUT. If N is the number of commutating 
capacitors in the network, ref. 2 shows that the following frequencies 
are present in the signals OUT and Y: (2N-1)w (2M1)w , (4N-1)u , 
(4N+1)a_ , é EC. 

In order to analyze the discontinous signal Y, a computer program 
was written. The object of the program was to calculate the components 


an and DO fied series OL the .torum 


a 
ew 
Y = 5 a 


it 18 


(a Cos mh 6 + Bb Sinema t) (2-9) 
1 Nn Nn n gi 


Nn 


ley) 


Figure 2.12 shows the program written in FORTRAN IV. Table I lists 
the component values a and DO of the output signal Y. 

It is seen from the table that the following components are dominant: 
Ist, 7th, 9th, 15th, 17th, 23rd, 25th, 3lst, 33rd, 39th, 41st, 47¢heeene 
etc. This is in agreement with theory given in ref. 2. 


An investigation was made to find out if the value of t the time 


0.” 
where Pi first goes to the value +1.0, had any effect on the output 
signal Y. The following phase delays were tried for the commutating 
wave-forms: 1/5, 27/5, 3n/5, and 4n/5. The results showed that phase- 


delaying the commutating square waves does not change the output signal. 


The series of equation (2-9) still has the same coefficients. 
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FORTRAN 


EME 


SUBROUTINE TORIT(FNT,N,M,A,B) 

DIMENSION A(100) ,B(100) ,FNT (10500) ,C(100) , PHI(100) 
WRITE (6,11) 

Bamore7x,'3'.8%,'A(J)', tore)’, 13X,°CCJ)' jet9xX, 


60 AN=N 


70 


75 


80 


12 


90 


100 


COEF=2 .0/ (2 .0*AN+1.0) 
CONST=3 . 141593 *COEF 
S1=SIN(CONST) 

C1=COS (CONST) 

P= 1.0 

S=0 ..0 

J=1 

FNTZ=FNT (1) 

U2=0.0 

U1=0.0 

T=2*N+ 1 

UO=FNT(I)+2 .0*P*U1-U2 
U2=U 1 

U1=U0 

I= I-l 

IF(I-1) 80,80,75 
A(J)=COEF *(FNTZ+P*U 1-U2) 
B(J)#COEF *S*U 1 

K = J-1l 

C(J) = SQRT((ACJ)) *(AGJ)) + (BYWI))* (BW))) 
PHI(J) = -ATAN(B(J)/A(J))*57.2957795 
WRITE(6,12) K,A(J),B(J),C(J), PHI(J) 
FORMAT (3X ,15,1P4E17.7) 
IF (J-(M+1)) 90,100,100 
Q=C 1*P-S1%*S 

S=C1*S+S 1*P 

Ee = a2 

J=J+1 

GO TO 70 

A(1)=A(€1)*0.5 

RETURN 

END 


Figure 2.12 Computer program for Fourier analysis 
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.0915914E-02 
.7432428E-02 
.439543 8E-04 
.7318883E-04 
.5144676E-04 
.0724834E-03 
.3770524E-03 
.2983000E-01 
.4599903E-04 
.9531114E-02 
.1894142E-03 
.7437004E-04 
.1180733E-04 
.8510741E-04 
.2539515E-04 
4.5283 84E-02 
.3075595E-04 
.28/9425E-02 
.3891549E-03 
.9694903E-04 
.3581335E-04 
.0385311E-04 
.0932467E-04 
.9760877E-02 
.1980087E-04 
.6020168E-02 
.6435747E-05 
.6750767E-05 
.2703399E-04 
.2502408E-04 
.3039294E-03 
.1548135E-02 
. 1624935E-05 
.0806191E-02 
.343 8673E-03 
.8415312E-04 
.2506500E-04 
.0879410E-04 
.7829986E-04 
. 198547 8E-03 
.8792193E-04 
.8365095E-03 
.8027469E-04 
.6715851E-04 
.9095364E-04 
.7992225E-04 
.4365583E-03 
.6978830E-02 
.2683252E-04 
.6212072E-02 
.0357 858E-04 
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0 
7068429E-01 


91272 Ee - 08 
.3683525E-03 
. 864587 8E-04 


2038115E-03 


947523 1E-05 


7617463E-02 


-4301525E-04 
.0481014E-02 


2319605E-03 


.9847650E-04 


1678812E-04 
2260592E-04 


-4881822E-03 


89047 10E-02 


.5192422E-04 


1087449E-02 
6604668E-04 
0924985E-04 
5659863E-04 
5742664E-04 


.7328471E-04 
-8151786E-02 
.3531574E-04 
.7231848E-02 
.4239037E-03 
.6518889E-04 
.1566226E-04 
. 1200222E-04 
.7783101E-04 
.6806425E-03 
.2581698E-04 
.378433 8E-03. 
.3319259E-04 
845423 9E-04 
-4045424E-04 
.9357740E-04 
-461943 8E-03 
-4380792E-02 
.8413758E-04 
770033 /9E-02 


1022575E-03 
9376331E-04 


~7252291E-04 


5532 831E-04 


.4710323E-04 
. 1989255E-02 


9340992E-04 
325598/7E-02- 
2116958E-03 


Table I Fourier coefficients of output signal Y 


IIIT. FREQUENCY RESPONSE OF FILTER 


An investigation of the phase difference between output and input 
signals was done for frequencies within one decade on either side of 
the notch gan ee wo: 

The frequency response was done for a notch frequency of one hertz 
and with a time constant of one second. 

The equivalent network of fig. 2.6 was used for the frequency re- 
Sponse study where the frequency of the input signal, x, was varied 
between 0.1Hz and 10.0 Hz. 

A Fourier analysis was done of theoutput signal y. The program of 
figure 2.12 was used for the analysis. The phase difference between 
the input signal and the first harmonic of the output signal was calcu- 
lated. The calculations were done for times when the phase of the 


input signal was zero. The phase, ¢ of the first harmonic of the 


17 
output signal is 


$= - tan (B,/A,) (3-1) 


Figure 3.1 shows the plot of phase versus frequency as obtained. 
The phase starts at zero degrees for low frequencies and goes to -90 
degrees for frequencies slightly less than Wo. jumps to +90 degrees for 
frequencies slightly above Ws and goes back to zero degrees for high 
frequencies. 


The values obtained in the computer simulation are collected in 


Table II. 
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Frequency response of filter - phase 


Figure 3.1 
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Frequency Phase 


1 ie Oh) 
OZ - 7.40 
0.25 - TT 
2 ee: “13.59 
0.4 -20.12 
OS -24.00 
0% -30.07 
Os -51.54 
0.2 <Fle23 
les 25 43.80 
2.0 15.14 
Pa 14.33 
4.0 4.38 
10.0 2.09 
20.0 0.28 


Table II Results of frequency response simulation - phase 
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IV. APPLICATIONS OF FILTER IN SERVOSYSTEMS 


Computer simulations were done with the filter used in conjunction 
with servosystems. The objective of the filter was to suppress the 
component of frequency in the output signal C. 

Example l 
A servosystem with the following transfer function was used in the 


first example: 


2 
c(s ) ee 


x1(s) ~ i 


7D 
s (0.00025 stl) GS = ro s+1) 
60) n 


n 
Figure 4.1 shows the system that has the open-loop transfer function 


of equation (4-1). 


X1i(s) + KHo: \ V Ci 
®; S+2W.S* We 0 AS+4 $ 


Figure 4.1 Servosystem of example l 
The following values were used: §€& = 0.01, wo, = anf 2500 = 15708 rad/sec. 
Figure 4.2 shows how the filter was inserted in the servosystem. Figure 


4.3 shows the Bode-plot for the system. 
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Figure 4.2 Compensated filter network 
The signal VN is the sum of the signals OUT and VM. Since the funda- 
mental component (frequency mca) of the signal OUT is 180 degrees out 
of phase with the signal VM, the sum signal VN will have no component 
of frequency @ . The computer simulation of the system of fig. 4.1 was 
carried for different values of the gain KMO. Figures 4.4 -~ 4.6 show 
the time responses of the system to a unit step input for different 
values of KMO. It is seen that the output of the system when the 
filter is inserted is much smoother than the output of the system 
without the filter inserted, and it does not contain any component of 
frequency @ . 

It is also seen that as the gain KMO is increased, the oscillations 
of the filtered output also increase in magnitude, but their amplitudes 
in the filtered output are still much smaller than their amplitudes in 
the unfiltered output signal. While the unfiltered output signal be- 
comes unstable, the filtered output is still stable. An analysis of 
the output signal C shows that the dominant frequency of the signal 
varies as a function of the gain KMO. The dominant frequency is found 
to be the frequency where the magnitude curve of the Bode-plot crosses 


the zero-decibel axis. 
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Example 2 - 4 


In examples 2 - 4 the same servosystem was used. Only the value of the 
gain KMO was changed. Figures 4.7 - 4.9 show the Bode-plots for the 
servosystem with and without the filter inserted. In example 2 the 
notch frequency is set at the phase cross-over (frequency = 15708 rad/sec). 
In example 3 the notch frequency of the filter is set at the magnitude 
cross-over of the Bode-plot. Finally in example 4 the notch frequency 
is set one decade below magnitude cross-over (frequency = 314 rad/sec). 

The computer simulations for the time responses for examples 2 - 4 
agreed with the expected results. 

In example 2 the unfiltered output is stable with a phase margin of 
about 13 degrees. The output when the filtered is in the network has 
less phase margin but is still stable. 

In example 3 the unfiltered output is stable with a phase margin 
of about 18 degrees. As seen from the Bode-plot for the system with 
the filter in the circuit the output is now unstable. 

In example 4 the unfiltered output is identical to that of example 
3. The notch frequency, however, is now one decade below the frequency 
w >» and the filter's frequency responses are now back to OQOdb and 0 
degrees for magnitude and phase plots respectively at the gain cross- 


over. The filter has thus no influence on the output signal. 


2] 


- | Wa << untilte ced cutpu 
| , 


| filtered outyut 











Figure 4.4 Time response to unit step input .KMO=8.636x10?+ 


Scale: X=0.0004 sec/in, Y=0.2 units/in. 
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Figure 4.5 Time response to unit step input .KMO=2.0x10* 


Scale: X=0.0004 sec/in, Y=0.2 units/in. 
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Figure 4.6 Time response to unit step input -KM=2/5x10"- 


Scale: X=0.0004 sec/in, Y=0.2 units/in. 
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V. APPROXIMATION OF FILTER TRANSFER FUNCTION 


As seen from fig. 3.1 the phase part of the frequency response is 
discontinous at the notch-frequency. 
A transfer function of the form 
Ba +o? 
1S) =e noe > (5-1) 
2 iz 
S Ge Za. s + 
n n 
is seen to have the same phase characteristic. By varying the value 
of the damping factor the curvature of the phase curve can be changed. 
Figure 5.1 shows the locations of the poles and zeros of the transfer 


function T(s) in the s-plane. The poles are shown for different values 


of the damping factor. 





Figure 5.1 Poles and zeros of T(s) 
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Figure 5.2 shows the Bode-plot of T(s) and of the commutating network 


filter as obtained in chapter III. 


35 


oseyd - (s)I Jo j0Td-epog zc ain3ty 





- 
ANA AE NIOS 08 GS, EN, PEM iN Le tT RT Lee ET ty rasan se 0% 8S De te TE EM ee et ers Seek ee te oe. 











OF re ae 


€°0=32 


36 


A9ITEF-NO 


= fe ee eet a oe 


te Ci age Te Fie rere AEE - Er hee ~ te mem GO nek tk ee ee OO EAS ee ee SE ye HS ee oe gy Fee Oe ee ee ee 


VI. CONCLUSION 


The commutated network filter has a highly selective bandpass 
filter characteristic where the center frequency of the bandpass 
frequency-response curve can be easily varied by changing the frequency 
of the commutating functions. 

The filter can be approximated by a linear second-order transfer 
function for design purposes. The transfer function has two zeros on 
or very close to the imaginary axis and two poles on a circle with 
radius w.. The location of the poles can be changed depending on the 


value of the damping factor. 
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APPENDIX A 


TIME RESPONSES OF SERVOSYSTEM 


A compensated servosystem with open-loop transfer function 


G(s) a= KMO 


= ———__ hn WL (A-1) 
xi (s) (ast1) (s*F2E w sto”) 


was used in example 1, chapter IV to get time responses to a unit step 
input for the system with and without the filter inserted. 

Additional simulations were carried out for different combinations 
of the natural frequency and the damping factor €. Figure A.1 shows 


which combinations of D and € were used. 





Figure A.l Values used for time responses 


Table III lists the values used for KMO, w and &, 


ers) 








Figure w /2n e xmo/10!+ X-scale Y-scale 
sec/in units/in 
Aap 2450 0.010 GS, ols 0.0004 Ov2 
Ae, 2550 0.01 8.463 0.0004 OZ 
A.4 2500 07005 8.636 0.0004 0.2 
A.5 2500 Beo15 8.636 0.0004 0.2 
A.6 2400 0.01 8.986 0.0004 0.2 
A./ 2600 GOL 8.290 0.0004 OsZ 
A.8 2500 0.0025 8.636 0.0004 OZ 
Pe? 2500 0.020 25.000 0.0004 G.3 
A.10 2500 0.020 8.636 0.0004 0.2 
A.ll 2250 0.010 9.499 0.0004 0.2 
A.12 2230 0.0025 9.499 0.0004 0.2 
Ay 13 2250 0.020 9.499 0.0004 OZ 
A.14 Z7 30 0.01 ETP 0.0004 Ore 
=. 15 2750 0.0025 fetie 0.0004 Oe2 
A.16 2750 C202 Lepee 0.0004 O52 


Table III Values used for computer-simulation 
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Figure A.2 Time response to unit step input 
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Figure A.3 Time response to unit step input 
/ 
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Figure A.4 Time response to unit step input 
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Figure A.5 Time response to unit step input 
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Figure A./ Time response to unit step input 
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Figure A.8 Time response to unit step input 
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Figure A.10 Time response to unit step input 
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Figure A.1l1 Time response to unit step input 
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Figure A.12 Time response to unit step input 
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Figure A.13. Time response to unit step input 





Figure A.14 Time response to unit step input 
4 





Figure A.15 Time response to unit step input 
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Figure A.16 Time response to unit step input 
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